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1 INTRODUCTION 



ABSTRACT 

The central engines of active galactic nuclei (AGN) contain cold, dense material as well 
as hot X-ray emitting gas. The standard paradigm for the engine geometry is a cold 
thin disc sandwiched between hot X-ray coronae. Strong support for this geometry in 
Seyferts comes from the study of fluorescent iron line profiles, although the evidence 
is not ubiquitously air tight. The thin disc model of line profiles in AGN and in X-ray 
binaries should be bench marked against other plausible possibilities. One proposed 
alternative is an engine consisting of dense clouds embedded in an optically thin, 
geometrically thick X-ray emitting engine. This model is further motivated by studies 
of geometrically thick engines such as advection dominated accretion flows (ADAFs) . 
Here we compute the reprocessed iron line profiles from dense clouds embedded in 
geometrically thick, optically thin X-ray emitting discs near a Schwarzchild black 
hole. We consider a range of cloud distributions and disc solutions, including ADAFs, 
pure radial infall, and bipolar outflows. We find that such models can reproduce line 
profiles similar to those from geometrically thin, optically thick discs and might help 
alleviate some of the problems encountered from the latter. 

Key words: accretion, accretion discs - line: profiles - line: formation - galaxies: 
active - galaxies: Seyfert - X-rays: galaxies 
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Active galactic nuclei (AGN) and some X-ray binaries are 
thought to be powered by accretion onto massive black holes 
(e.g. Rees 1984; Tanaka 2000). Observations suggest that 
cold matter in the accretion engine reprocesses the primary 
X-ray continuum (e.g. George and Fabian 1991). In partic- 
ular, a fluorescent iron Ka emission line profile is imprinted 
on the spectrum near 6.4 keV. The iron line profile is shaped 
by Doppler and gravitational frequency shifts (Fabian et al. 
1989), and hence contains information about geometry and 
dynamics in the innermost regions of the accretion engine 
(see Fabian et al. 2000 for a comprehensive review). 

Most modelling of iron line profiles has been carried out 
within the paradigm of geometrically thin, optically thick 
discs. The initial flat disc model of Fabian et al. (1989) was 
generalised to Kerr black holes by Laor (1991). More re- 
cently, Pariev and Bromley (1998) have considered effects of 
a small disc thickness and turbulence. 

Thin disc line profiles successfully reproduce many fea- 
tures in the experimental data (e.g. Tanaka et al. 1995; 
Nandra et al. 1997). However there are some observations 
which remain difficult to understand within this paradigm. 
For example, a preference for low inclinations, even amongst 



Seyfert 2s (Turner et al. 1998), disagreements between the 
predicted inclination from iron line modelling and the pre- 
dicted inclination from other phenomena (Nishiura et al. 
1998; Sulentic et al. 1998a; Wang et al. 1999), and problems 
with the observed reprocessed fraction being too low or not 
correlated with line width (e.g. Lee et al. 2000; Chiang et 
al. 2000; Weaver 2000; Sulentic et al. 1998b). 

Such puzzles may be resolved within the paradigm of 
thin discs, for instance by introducing additional emission 
components from an obscuring torus (Weaver and Reynolds 
1998) or invoking flared or warped discs in which the inner 
and outer regions are not aligned (Blackman 1999; Hartnoll 
and Blackman 2000). In any case, alternative plausible ge- 
ometries still merit exploration to provide proper diagnostics 
for comparison. 

A particularly interesting alternative extends from work 
on the survival of cold clouds capable of reprocessing X-ray 
flux in the AGN engines (e.g. Ferland and Rees 1988; Guil- 
bert and Rees 1988; Celotti, Fabian and Rees 1992). Phys- 
ical constraints on the size of the dense clouds have been 
studied by Kuncic et al. (1996). They found a sizeable pa- 
rameter regime in which cool, dense, small clouds can coexist 
with a hot medium in a dynamical equilibrium. Though X- 
ray continuum features of the embedded cloud model have 
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been studied by several authors (Sivron and Tsuruta 1993; 
Bond and Matsuoka 1993; Nandra and George 1994; Collin- 
Souffrin et al. 1996), detailed calculations of reprocessed line 
profiles have not been carried out. Karas et al. (2000) have 
studied line emission from reprocessing by a nearly spherical 
shell of clouds at a specific radius that results from a stan- 
dard thin disc being disrupted due to instabilities. They give 
resulting line profiles from illumination by a central source. 
This is a different scenario from what we will consider here. 
Other related calculations of note include Cao and Zhang 
(1991), who calculated line profiles for radially infalling gas. 

In this paper we calculate the (iron) line profiles that 
result from various distributions of small, cold clouds em- 
bedded in an optically thin, geometrically thick accretion 
disc or corona. As well as illumination from a central source, 
we include an illumination by an 'X-ray bath' of radiation; 
each cloud is taken to be small enough to be illuminated 
isotropically at their positional location by this bath com- 
ponent. This two component illumination is motivated from 
considering geometrical thick models such as, for example, 
advection dominated accretion flows (ADAFs) (e.g. Rees et 
al. 1982; Narayan and Yi 1994; see Narayan et al. 1998 for 
a review), which can be interpreted to have a synchrotron 
+ Bremsstrahlung centrally dominated illumination (Ma- 
hadevan 1997) and a Bremsstrahlung X-ray bath. Depend- 
ing on the accretion rate of the disc, the synchrotron or 
Bremsstrahlung will dominate. ADAFs predict luminosities 
significantly less than the Eddington value Le, and may just 
be applicable to Seyferts, with typical luminosities L_e/100. 

We will use several flow solutions to prescribe a velocity 
field for the clouds; standard full ADAF solutions, pure ra- 
dial infall, and bipolar outflows. The profiles calculated may 
also therefore be relevant for the dynamics of convection 
dominated accretion flows (CDAFs) (Quataert and Gruzi- 
nov 2000; Narayan et al. 2000), Advection Dominated Inflow 
Outflow Models (ADIOS) (Blandford and Begelman 1999) 
or other thick disc models. The microphysics of the direct 
X-ray line emission for ADAFs have been studied in non- 
relativistic flows (e.g. Narayan and Raymond 1999; Perna et 
al. 2000), but not the presence of cold reprocessing matter, 
nor the associated relativistic iron flourescence line profiles. 
Like Karas et al. (2000), we take the clouds to be optically 
thick, and allow the clouds to shadow emission from other 
clouds. 

Some profiles for relativistic outflows have been calcu- 
lated by Wang et al. (2000). These outflows are different 
from those considered here in that they are pencil beams 
along the spin axis and not radial outflow in a geometrically 
thick disc. The spherical geometry, and covering by clouds, 
can avoid the large blue shifted component that has not been 
observed in iron lines and that is characteristic of outflows. 

It has been suggested (e.g. Nandra 2000) that a geo- 
metrically thick distribution of clouds would result in more 
blue-shifted emission than is observed. The data from a va- 
riety of sources are also not inconsistent with a somewhat 
larger blue wing (Nandra et al. 1997). But we find that be- 
cause models of geometrically thick discs can generally have 
smaller bulk velocities in the central regions than the flat 
disc keplerian velocities, smaller Doppler blue shifts result, 
and gravitational redshift becomes the dominant spectral 
feature. For these reasons also, the line profiles for thick 
discs will in general be less inclination dependent than those 



from thin discs. This will be discussed in more detail below. 
It may help explain problems encountered with inclinations 
in thin disc modelling. When the velocities are faster, the 
profiles are very similar to the standard profiles obtained 
from thin disc modelling. 

In section 2 we describe the disc model of dense clouds 
in geometrically thick discs. Section 3 contains the calcula- 
tion of iron line profiles from this model. Section 4 presents 
the resulting line profiles and a discussion of the main fea- 
tures. Section 5 is the conclusion. 



2 THE EMBEDDED CLOUD MODEL: DENSE 
CLOUDS AND THICK DISCS 

The model consists of small optically thick clouds embed- 
ded in, and hence following, the fluid flow in an other- 
wise optically thin disc. Although the clouds may have a 
short lifetime, the continuous creation and destruction of 
clouds allows the definition of an effective number density 
of clouds at each point, n c (r). We consider these clouds to 
be both bathed in an locally isotropic radiative flux (e.g. 
Bremsstrahlung) and to be illuminated by a central source 
(e.g. synchrotron dominated). The clouds are assumed to be 
small compared to disc length scales (see Celotti et al. 1992 
and Kuncic et al. 1996 for partial justification). 

Calculating the profile then involves an optical depth 
calculation; we integrate the number density times the cross 
section along the photon trajectories to obtain a probabil- 
ity of obscuration. Contributions are integrated across the 
whole disc to obtain a profile, as is done for the usual thin 
disc, except that now we are integrating over three dimen- 
sions. 

We consider several different solutions to the accre- 
tion flow involving optically thin, geometrically thick en- 
gines: self-similar ADAFs (Narayan and Yi 1995; Perna et 
al. 2000), bipolar outflows (Xu and Chen 1997; Perna et al. 
2000), and radial infall (appendix B of Narayan and Yi 1995, 
similar to Bondi 1952). 

2.1 self-similar ADAF 

Self-similar accretion solutions are described in Narayan 
and Yi (1995). The approximate analytic solutions, given in 
Perna et al. (2000), for the density (p), velocities (v r , vg, v^), 
and sound speed (c 3 ) are 



P = p r-^[c%e)]- 2 - 2 \ 

3a _ . . . 2 
• - r\lK\'r) sin ( 



5 + 2e 



Vg = 0, 



V,f, = 



2e 



,5 + 2e 
c s = rQ, K {r) 



rQ.K{r) sm9, 



(1) 
(2) 
(3) 

(4) 
(5) 



-5 5 + 2e, 
where Q.k(t) is the Keplerian angular velocity, Q.k(j) = 

_ 3 

r 2 in our units and we are using standard spherical polar 
coordinates, (r,d,4>). 

There are two free parameters, the Sakura-Sunyaev vis- 
cosity parameter, a, and a thermodynamic parameter, 
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5/3-7 



7- 



(6) 



with 7 the ratio of specific heats. We will consider a — 
0.35 and e = 10,0.1. Solutions are taken to be in the fully 
advection dominated case. Note that: 

(i) Lower e gives lower angular speeds. 

(ii) For e — 10 most of the matter is close to the equatorial 
plane, whilst for e = 0.1 the matter is essentially uniformly 
distributed in 9. 

(iii) For larger values of a, such the a — 0.35 we consider, 
v r is very similar to v^. For smaller a, the radial motion 
becomes less important. 

(iv) c s (6) is smaller on the equatorial plane and is smaller 
for larger e. This will mean that, at fixed radius, the clouds 
are larger on the equatorial plane and larger in high e cases. 

2.2 Radial infall 

The self-similar forms of the previous section admit a radi- 
ally infalling solution (Narayan and Yi 1995, it is essentially 
Bondi accretion with viscosity.) We take the radial infall to 
be at the Keplerian velocity, rfljf(r), which implies that the 
sound speed is c s = r *L ■ This solution corresponds to 
a small e, and so represents a geometrically thick, optically 
thin, advection dominated flow. The density is taken to be 

_ 3 

constant in 6, so p = por 5 . 



2.3 Bipolar outflow 

Outflow solutions have been considered by several authors 
(e.g. Xu and Chen 1997; Quartaert and Narayan 1999; 
Blandford and Begelman 1999). Perna et al. (2000) give ap- 
proximate analytic expressions. We use their solutions with 
density independent of angle and with winds of p — | , where 
p is defined so that M oc r v . The equations then become 



P = Por 2 , 

3q 



rfiK(r)[cos20 + cos 2 0], 



Vg — 



W0 = 



5 + 2e 

3a rf2isr(r) 
5 + 2e 4 

2e 



sin 26», 



5 + 2e 
rVL K {r) 



rQ.K(r) sin( 



v / 2-ecos 2 6'. 



(7) 
(8) 

(9) 

(10) 
(11) 



3 CALCULATION OF THE LINE PROFILE 
3.1 Total flux 

The total flux observed at frequency v is given by 



F(y) 



5{v-v (r))f(r,v)d r, 



(12) 



where fo(r) is the observed frequency of flux emitted from 
r, calculated below, and 7(r, v) is given by 



/M = 



/ «t(r,c(r))V(r), 



(13) 



where y£-\ is the relativistic correction to the flux, with 
the frequency of emission, p e — 6.4keV for the iron line. 
foutijc, c(r)) is the flux reprocessed and re-emitted by clouds 
at r in the direction c(r), defined such that the photons reach 
the observer (calculated below). Finally, ^(r) is the visibility 
function which gives the probability of flux from r reaching 
the observer without being intercepted by an optically thick 
cloud. It is given by 



V(r) = max(l-P(r),0), 

with the probability of interception 



P(ro) 



n c (r(A))A(r(A) 



dv N 
lx' 



dr 
d\ 



dX, 



(14) 



(15) 



where r(a) = ro, the position of emission, and r(6) = r bs, 
the position of the observer, and A parametrises the path 
of the photon. n c (r) is the number density of clouds in the 
disc and A(r, 4f) is the cross sectional area of the clouds 
at r, as seen by a photon travelling in direction 4v. These 
quantities are discussed below. The visibility function must 
also include the possibility of photons falling into the event 
horizon. 

The dimensions of the disc are taken to be 6 < r < 300, 
< 6 < n and < <f> < 2n, although there will not necessar- 
ily be a significant number density of clouds throughout this 
whole region (because the number density of particles in the 
disc itself gets smaller for larger r and, depending on the 
thermodynamic parameter e, gets smaller further form the 
equatorial plane). Also, the volume element has a general 
relativistic correction, so 



dr 



V^ 



(16) 



3.2 Photon paths 

We use paths in the Schwarzchild metric, approximated by 
a first order perturbation. We take units with c = G = 1 
and take the mass of the central object to be M — 1. Note 
that in this section r is the position of the photon travelling 
from the cloud to the observer, and ro is the position of the 
emitting cloud in the disc. We write ro =| ro | and r =\ r |. 
First order perturbation of the Schwarzchild null 
geodesic equations gives the photon path to be 



r($) 



sin $ 1 + C cos $ + cos 2 $ 



2R 



AR? 



Ot 



i? 3 " 



(17) 



Where R and C are constants, r(<fr) is the radial distance 
from the origin at the central object and $ is a polar angle on 
the plane of the trajectory. The trajectory plane is defined 
by the origin, the direction to the observer, n, and the initial 
position vector ro. 

The boundary conditions are that the photon starts at 
(ro, $o) an d ends at (oo, $ 6 S ), written in (r, <3>) coordinates 
on the trajectory plane. The boundary conditions are sim- 
plified if we choose the axis on the trajectory plane that 
defines $ to be in the direction towards the observer, n, so 
that $o6 S = 0. The constants are then given by 



C=-2, 



_± ro sin $o ro 
ti = ■ ± 



.,sin 2 $ H (l-cos$ ) 2 . 

1 V r 



(18) 



(19) 
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Note that the minus sign case results in R negative, which is 
not physical, so we take the positive solution. Selecting this 
single path is effectively not considering orbits that loop 
around the central hole one or more times, consistent with 
the perturbative approach. 

With this choice of axis, "J>o is given by 



cos«l>o 



r n 

ro 



(20) 



In principle, there could be a problem inverting this expres- 
sion due to ambiguities in the definition of the inverse cosine. 
However, these are avoided if we take the other axis on the 
trajectory plane to be r °~ (r °' n)n . This vector, orthogonal to 
n, is such that r has positive component in this direction. 
Therefore, we may take $o S [0, n] with no ambiguities and 
the trajectory is given by 



[•($) = r($) 



cos <J?n + sin <E> 



r - (r ■ n)n 



'■o 



This gives the probability of interception as 
P(ro) 



dr. 



n(r(*))A(r($),-— ) 



dr 



and the angle of emission is given by 
rfr 



c(r ) 



d$> 



($ = $o) 



(21) 



d$, (22) 



(23) 



The minus sign in the last two equations results because the 
photon travels in the direction of decreasing $. 



3.3 Frequency shifts 

The observed frequency is given in terms of the emitted 
frequency by 



i/ = v t 



(P^obs 



(24) 



where p^ is the photon 4-momentum and £/ M is the 4- velocity 
of the emitter or observer. We employ general relativity in 
the Schwarzchild metric, g^. Write p M = (p',p) M , and £/ M = 
d7(l,u) // - The time component of photon 4-momentum in 
the Schwarzchild metric is 
E 



P = 



1 



(25) 



Where E is constant along the photon trajectory. So requir- 
ing the photon 4-momentum to be null, g^vP^p" = 0, gives 
P 2 — QijP l p' — rrr to be the appropriate normalisation for 
the spatial components. In the Schwarzchild metric, we have 
dt _ 1 

dr rr 



(26) 



By definition, 
Pv-U 1 * = g^p^U" 



(27) 



where the observer has u (, s = and r (, s — > oo. The emitter 
has p em oc c and r em — rg. Putting all this into (|24J) gives 



vo (ro) 





-V 1 -*- 


_^_ M 2 (ro) _ M 2 (ro) 
r 


1- 


c r (r )u r (r ) 

1-JL 


- ce(ro)ue(ro) — c^(ro)it^,(ro) 



(28) 



Here u(ro) is the velocity of the clouds at ro. We take the 
velocity of clouds to trace the velocity of the fluid in the 
disc. This is given by one of the disc motions described in 
the previous section. 



3.4 Distribution and sizes of clouds 

For simplicity, we will assume that all of the clouds are 
spherical. The cross section is then independent of a pho- 
ton's direction of approach. If we consider clouds with the 
same total number of particles, this allows us to to write 



A(r) oc Va(r) oc 



= W 



(29) 



where n PjC is the number density of particles in the cloud. 

Thermal equilibrium gives a relation between the num- 
ber density of particles in the cloud and the number density 
of particles in the disc, n, 



n(r)T(r) 

Tc(r) 



(30) 



where T and T c are the temperatures of protons in the disc 
and the clouds, respectively. We take the temperature of the 
clouds to be the same everywhere - around T c — 10 5 K (Kun- 
cic et al, 1996), although the actual value will be absorbed 
into an overall constant - and find the temperature of the 
disc from the sound speed, c s , as 



T(r) = 



m p c 2 s (r) 



(31) 



where k is Boltzmann's constant and m p is the mass of the 
proton. This gives us the cross section 



A(r) oc 



kT c 



m p n(r)c|(r) 



p(r)c!(r) 



(32) 



We take p(r) and c^(r) from the solutions of the previous 
section. The expression should be normalised by the cross 
section of the clouds being around 10 10 cm 2 at the inner- 
most radius (Kuncic et al., 1996), however in practice we 
can absorb this constant into the constant F, defined below. 
The various assumptions above about the cloud struc- 
ture (constant number of particles, constant temperature) 
result in clouds with size growing strongly with r, as r'. 
Combined with the radial dependence of the flux emission 
(discussed in more detail below), this results in profiles dom- 
inated by the outer disc regions. We also want to consider 
more general radial dependencies, where profiles are dom- 
inated by the central regions. We therefore add the term 
r~ g , to the cloud cross section, where generally g is a small 
positive number. This has a somewhat analogous role to the 
r~ q term in the standard thin disc models (Fabian et al. 
1989), although has a physically different origin. The final 
expression is then 



A(r) oc 



p(r)c!(r) 



(33) 



We take the number density of clouds to trace the num- 
ber density of protons in the disc, hence 



n c (r) = Tp(r). 



(34) 
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The constant of proportionality, V, will determine the degree 
of shadowing by optically thick clouds. The probability of 
obscuration, P(r), with maximal value 1, is proportional to 

r. 



3.5 Emitted flux 

We take the emitted flux to be proportional to the flux illu- 
minating the clouds, /;„, the number density of the clouds, 
and the cross section of the clouds. 

We take two contributions to the flux illuminating the 
clouds, a central illumination (e.g. due primarily to syn- 
chrotron radiation in the innermost regions of the flow for 
ADAFs Mahadevan 1997) an X-ray bath throughout the disc 
(due to Bremsstrahlung radiation). The central illumination 
falls of as 4j and the Bremsstrahlung is proportional to the 
square of the number density of particles in the disc. So the 
emitted flux is given by 



fout(r, c(r)) oc n c (r)A(r)/in(r) 

., imaxf- f • c(r),0) 
oc pA(r)3 



V p (r) + Bp 2 



(35) 



Here A and B are constants determining the relative im- 
portance of central radiation and the X-ray bath. They 
will depend on the mass accretion rate of the system (Ma- 
hadevan 1997) and also the geometry of the clouds. The 
max(— f -c(r), 0) term contains the information that only the 
side of the cloud facing the centre of the disc is illuminated 
by the central flux. Hence, the corresponding reprocessed 
flux is only visible for clouds on the far side of the hole with 
respect to the observer, and the amount of flux seen by the 
observer decreases as the illuminated side of the clouds be- 
comes more edge on to the observer. V p (r) is the primary 
visibility function, analogous to the visibility function de- 
scribed earlier, giving the probability of photons emitted 
from the central regions reaching the cloud without being 
obscured by intermediate clouds. Fortunately, the photon 
paths from the central illumination are radial and so the in- 
tegration of equation Il5| becomes analytically tractable, us- 
ing appropriate expressions for the cross sections and cloud 
number densities. 

To give some visualisation of the processes involved, we 
will consider five cross sections of the disc at different heights 
above the equatorial plane, as shown in Fig.l. For each cross 
section, we compute a frequency map and a visibility map 
for both a centrally illuminated (synchrotron dominated) 
case and an X-ray bath (Bremsstrahlung dominated) case. 
The disc visualised is a self-similar ADAF disc with e = 0.1 
and a = 0.35, viewed at 60 degrees. The results are shown 
in Fig. 2, in which the observer is situated to the right. 

Some thought reveals that the plots of Fig. 2 make sense 
physically, and the qualitative features deserve some further 
comments. The redshift map of the equatorial is different 
in two respects from the standard redshift maps of thin 
discs. Firstly, the regions of largest Doppler shift are ro- 
tated slightly with respect to the observer. This is because 
the radial motion is of similar magnitude to the rotational 
motion. Secondly, gravitational redshift is much more domi- 
nant. This is because the velocities are significantly less than 
the Keplerian velocities of a thin disc (see comments in §2.1). 




Figure 1: The cross sections of the disc used in Fig. 2. 



As we would expect, the redshifts become less pronounced 
as we move away from the equatorial plane. 

The visibility in the centrally illuminated case is a su- 
perposition of the visibility in the X-ray bath case with the 
primary visibility. The primary visibility implies that the ob- 
server can only see clouds on the far side of the hole (hence 
there is more obscuration on the top part of the disc than 
the bottom part) and that the more edge-on clouds appear 
less illuminated. 

The visibility in the X-ray bath case is determined by 
the size and number density of clouds as well as the the 
length of the of the trajectory through the clouds. This re- 
sults in the central regions being the most obscured because 
the disc under consideration has g — 1.5, which means the 
cloud number density times cross section is largest in the 
central regions. Also, the bottom regions are more obscured 
than the top regions because they are further from the ob- 
server. 



4 LINE PROFILES 

Discs were considered (i) at two or three different inclina- 
tions (15,30 and 60 degrees, or some subset thereof), (ii) 
once with shadowing (r finite) and once without (r = 0), 
and (in) once with the profile dominated by the outer re- 
gions (g — 0) and once with the profile dominated by inner 
regions (g — 1.5). All discs have a = 0.35. 



4.1 Self-similar ADAF profiles 

We consider three ADAF discs. The first has e = 0.1, hence 
the density dependence is essentially spherically symmetric. 
We take this disc to be dominated by central illumination 
(i.e. A — 10, B = 1 in equation pq). The second disc also 
has e = 0.1 but is dominated by the X-ray bath (i.e. A — 
1,B = 100). The third has e = 10, hence the density is 
concentrated within around 30 degrees of 6 on each side of 
the equatorial plane. We take this disc to be dominated by 
central illumination. 

A striking feature of the profiles from the two e = 0.1 
ADAF discs considered, Figs. 3 and 4, is the lack of signif- 
icant dependence on the inclination angle (compare to the 
profiles from flat thin discs in, e.g. Fabian et al. 1989, which 
show a strong inclination dependence). This is because the 
velocities given by equations Q-M have magnitudes smaller 
than the Keplerian velocities of thin discs and so extreme 
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Figure 2: From top to bottom, cross sections of a disc at heights of 0.7r ou t, 0.25r ou t, 0, — 0.25r ou t and — 0.7r ou t vertically relative to 

the equatorial plane, drawn to scale. The left column shows redshifts (darker shows more redshifted), the central column shows 

visibility in the centrally illuminated (synchrotron dominated) case and the right column shows visibility in an X-ray bath 

(Bremsstrahlung dominated) case. White regions are totally unshadowed and black regions are totally shadowed. The observer is to the 

right and the disc is rotating counter clockwise. 
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Figure 3: Profiles from an ADAF disc with e = 0.1. Central 
illumination dominated (A = 10, B = 1). Left to right: non- 
shadowed and shadowed. Top to bottom: g = (outer region 
dominated) and g = 1.5 (inner region dominated). Each plot 
shows profiles for inclinations of 15 and 60 degrees, plotted to 
scale (i.e. not separately normalised). The profile with the higher 
peak is at 15 degrees. 



Figure 4: Profiles from an ADAF disc with e = 0.1. X-ray bath 
dominated (A = 1, B = 100). Left to right: non-shadowed and 
shadowed. Top to bottom: 3 = and g = 1.5. Each plot shows 
profiles for inclinations of 15 and 60 degrees, plotted to scale (i.e. 
not separately normalised). The profile with the higher blue cutoff 
is at 60 degrees. 



blue Doppler shifts do not appear. There is more blue-shifted 
flux for the higher inclination discs, which is why the peak 
very near the rest frequency is lower in the higher inclination 
discs. 

The discs with g = 1.5 have profiles dominated by flux 
from the inner regions, which explains the greater presence 
of gravitationally redshifted flux. This is most notable in the 
X-ray bath case of Fig. 4 because here there is a contribution 
from Doppler redshifted flux which is not so visible in the 
centrally illuminated case. See Fig. 2 for how the shadowings 
and redshifts interact. The fact that v r ~ vq in the discs of 
Figs. 3 & 4 is important because it rotates the most Doppler 
redshifted region into the less visible sector of the disc for 
the centrally illuminated case, Fig. 3. 

The most notable feature of the profiles from the e = 10 
ADAF disc, in Fig. 5, is that cases with g — are of virtually 
identical form as the standard thin disc profiles, and show 
the same inclination dependence as the thin discs! This is 
not too surprising because the velocities of equations H-ph give 
a magnitude very close the Keplerian magnitude and the 
angular motion is more important than the radial motion. 

The profiles with g — 1.5 are dominated by emission 
from the central regions, which results in the large amounts 
of red and blue shifted flux due to Doppler and gravitational 
effects. 
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Figure 5: Profiles from an ADAF disc with e = 10. Central 
illumination dominated (A = 10, B = 1). Left to right: non- 
shadowed and shadowed. Top to bottom: g = and g = 1.5. 
Each plot shows profiles for inclinations of 15, 30 and 60 degrees, 
plotted to scale (i.e. not separately normalised). The profiles with 
higher blue cutoffs are at higher inclinations. 
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Figure 6: Profiles from a radially infalling disc. Central illu- 
mination dominated (A = 10, B = 1). Left to right: non- 
shadowed+shadowed (to scale) and shadowed (scaled up). Top 
to bottom: g = and g = 1.5. 



Figure 7: Profiles from a bipolar outflow disc. Central illu- 
mination dominated (A = 10, B = 1). Left to right: non- 
shadowed+shadowed (to scale) and shadowed (scaled up). Top 
to bottom: g = and g = 1.5. Inclination of 30 degrees. 



4.2 Radial infall profiles 

We consider a radially infalling disc dominated by central il- 
lumination. The profile is independent of inclination because 
of the spherical symmetry, so only one inclination angle is 
considered. The profiles are shown in Fig. 6. To show ex- 
plicitly the effect of shadowing, the non-shadowed profiles 
are plotted to scale with the shadowed profiles in the left 
column. 

The profiles show a significant amount of blueshifted 
flux and little redshifted flux. This is the immediate con- 
sequence of the central illumination which means the blobs 
falling inwards towards the observer are visible, whilst those 
falling inwards away form the observer are illuminated on 
the side obscured from the observer. In the case where the 
profile is dominated by contributions from the central re- 
gions (g — 1.5), we see a gravitational redshift and more 
blueshift from the fast inner flow. 



4.3 Outflow profiles 

We consider a disc with the bipolar outflow solutions of 
equations U\ to Pill The disc has e = 0.1 and is dominated by 
central illumination. 

The profiles are very similar to those generated by the 
ADAF solutions without outflows. Here we show the 30 de- 
gree inclination profiles, but the profiles for different incli- 
nations are very similar, as for the ADAF cases. The double 
peak for the outer region dominated case is seen in the top 
row of Fig. 7. The redshift dominates the flux in both cases. 
This is because of the gravitational redshift, and because the 
use of central illumination means that the observer primar- 



ily sees emission from the clouds moving away from his/her 
line of sight. This is opposite to the radial infall case. 



5 DISCUSSION AND CONCLUSION 

We have considered a model in which cold, dense clouds are 
embedded in an otherwise optically thin, but geometrically 
thick accretion flow. Assuming these clouds reprocess pri- 
mary X-ray emission from a central source and an X-ray 
bath, we have computed the expected fluorescent iron line 
profiles. This is a different scenario to the standard calcula- 
tion of line profiles from thin, flat discs. 

We considered various solutions for the dynamics of the 
disc: ADAFs, pure radial infall models, and bipolar outflows. 
We allow these solutions to specify the velocity, number den- 
sity and size of the clouds at different points of the engine. 
In calculating the line profiles, a new effect compared to the 
thin disc models is the obscuration of flux by the clouds 
themselves, which are optically thick. 

Many of the resulting line profiles are similar to those 
from thin disc models, suggesting that thin discs do not pro- 
vide a unique model that is able to reproduce observations 
of iron line profiles. This suggests that cloud models warrant 
further study. Previous work has focused primarily on the 
X-ray continuum (Sivron and Tsuruta 1993; Bond and Mat- 
suoka 1993; Nandra and George 1994; Collin-Souffrin et al. 
1996) but there has been little modelling of the relativistic 
iron fluorescent line for the cloud model. Karas et al. 2000 
considered non-relativistic profiles with large clouds at one 
radius. 

A further incentive to consider thick disks with embed- 
ded clouds is that they may help resolve problems arising in 
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thin disc modelling which relate to the inclination of discs 
and the reprocessed fraction. This includes the puzzling pref- 
erence for low inclinations for Seyfert 2s (Turner et al. 1998), 
some incongruences between inferred inclinations from iron 
line modelling and those determined from other methods 
(Nishiura et al. 1998; Sulentic et al. 1998a; Wang et al. 
1999), and the observed reprocessed fraction appearing to 
be too low and/or uncorrelated with line width (e.g. Lee et 
al. 2000; Chiang et al. 2000; Weaver 2000; Sulentic et al. 
1998b). In this regard, we find that many of our embedded 
cloud, thick disc profiles (modelled with ADAF type flows) 
are relatively insensitive to the inclination angle because a 
significant pressure support means that the bulk velocities 
are smaller than the Keplerian values. 

This latter feature also means that many of our pro- 
files have sharp blue cut-offs because they are dominated 
by the gravitational redshift. Whether the sharp blue cut- 
off remains true in the Kerr case, where the disc extends 
closer to central hole and hence can acquire greater veloci- 
ties should be investigated. Note that this will also depend 
fundamentally on the innermost radius at which the repro- 
cessing clouds can survive. A thick geometry can also allow a 
wider range of reprocessed fractions, depending on the cloud 
distribution and amount of obscuration. 

To summarise, the embedded cloud+ thick disc model 
allows a wider range of disc flow solutions to show repro- 
cessed X-rays lines, and highlights how an alternative ge- 
ometry can produce profiles similar to those from thin discs. 
The model can also lead to profiles with different properties 
to those of the thin discs, such as less sensitivity to incli- 
nation angle. As such, we suggest the cloud model deserves 
further investigation. This would include more consideration 
of the mechanisms of cloud formation and sizes, generalisa- 
tion to the Kerr metric, and application to specific sources. 
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